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Abstract 

Bernoulli numbers are usually expressed in terms of their lower in- 
dex numbers (recursive). This paper gives explicit formulas for Bernoulli 
numbers of even index. The formulas contain a remarkable sequence of 
determinants. The value of these determinants for variable dimension is 
investigated. 
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1 Introduction. 

In a previous short paper (ref. arXiv: math. NT/0503160 vl 8 Mar 2005) the au- 
thor presented an explicit formula for Bernoulli numbers of even index based on 
a sequence of determinants. This particular formula was an accidental byprod- 
uct of a different mathematical investigation. The proof via this route was fairly 
tedious and it was decided to just publish the formula and see wether there was 
any interest in it at all. 

Mr Robin Chapman at the University of Exeter in the UK promptly re- 
sponded and provided a short and elegant proof using geometrical transform 
techniques (generating functions). Also a request was received from the Math- 
ematical Society in Korea for providing a proof for their Newsletter of JMS. 
Further requests were received from China. 



As a result this new version of the earlier short paper provides three similar 
explicit expressions for Bernoulli numbers. A first one is derived based on 
the sequence of determinants mentioned earlier. The proof of the second and 
third one starts from the first one and uses transform techniques similar to Mr. 
R.Chapman's. Finally the value of the sequence of determinants is investigated. 
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2 Non-recursive expression for even-index Bernoulli 
numbers. 

First a sequence of determinants Dk (with k= 1, 2, 3, ...) is denned as follows: 
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For convenience we also define 

Co = 1 (2.1) 
Let B 2p (with p= 1, 2, 3, ...) be the even Bernoulli number of index 2p. 
B 2p is then given by either of the following formulas: 
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3 Relation between the determinants . 

The determinants Dk in (2.1) are in explicit form. Nevertheless a recursive 
relation among the D^'s may be given as is obvious from their definition. This 
relation may be written successively as: 

D 1 = i D 

D2 = 51 D\ - Do 

D 3 = i D 2 - i D 1 + i Do 
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or in general: 



Equivalently: 



§7^ = ° (3 ' 2) 



4 Proof of formula (ED- 



We start from the well-known recursive relation among even-index Bernoulli 
numbers: 



p 
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^ I 2k - 1 y 2fc 2 2p + 1 
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The easiest way to proceed is to take a specific case, e.g. start with p=l, 2, 3: 



2\^ = 1_1 
1 J 2 ~ 2 3 



4\^2 /4\B4_1_1 
1 J 2 + I 3 J 4 ~ 2 5 



6\-B2 /6\S ± /6\ : B6 = l l 
1 2 + 3 4 5 6 2 7 



After some manipulation and in matrix form: 
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The determinant of this system equals 1 . 
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Solving for 5f and keeping in mind the definition of the determinants 
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It is obvious that this procedure may be continued for ever increasing values of 
p. This proves Q2.2p. 

5 Proof of formulas (E3D and (1211) . 

Notice that formulas i|2.2|) . I|2.3|) and (|2.4|l also apply to the case p=0, i.e. B =l. 
In terms of geometric transforms, parts of formulas l|2.2|l and (|2.3|) come close 
to a convolution. 

By expressing the transforms of both (^jyr an d -^p as series in even powers 
of the complex variable z we achieve convolutions. 

Let: 
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Convolving -D(^) with F(z) 
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G{z)=D(z)F{z) = Y,Y. 
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Since (see (J22J>): 

Sf=o ^~(2i+i)\~ l = ^ or k=l,2... and 1 for k=0, the result is the unit impulse 
i.e. G(z)=l. ' 
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Therefore: 



D(z) = 4- (5.2) 



Looking now at formula ^2.2(1 . the geometric transform of 

J(Z)C ° SZ + J(Z) = D (*) cos 2 H = Zg8 ! f ! (5-3) 
2 w V2/ 2sin(f) V 7 

By ()2.3|1 the transform of the same expression 

"W = "W + 2* g (20! ° qUalS: 
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z (4 sin 2 (f ) cos (§) + 4 cos 3 (f) - 3 cos (f )) zcos (§) 
= 2sin(f) = 2sin(f) 

i.e. Q5.3J1 and l|5.4|l are identical. This proves formula (|2.3(l . 

We proceed in a similar way for formula (|2.4|) : 

Let: W(p) = 

and: W^) = E~ o WV 2p 
We already know: 

, % zcos (§ ) 
2sm(f) 

Formula |HI> claims: D p = - (2 2 p - 2) W p 
The geometric transform of (|5.5|l equals: 

1 /-2zcosz 2zcos(|)\ z 

2 I sinz sin (|) y sin 

This proves l|2.4ll . 
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6 Value of the determinants D^. 

From H5.2fl : 



D{z) = ^- = ^D P ^ (6.1) 
sin z 

p=0 



»v=-i^ = -i^- (6-2) 

where the integral is taken along a sufficiently small circle around the origin 
z=0. 

We now use a well-known result: 
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For z| (7r and k=l,2,.... , we are allowed to expand the last term in a power 
series. By dividing further by z 2p we obtain: 
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The value of l|tj.2fl equals the residue of . 

For p=0 the residue equals 1, i.e. D — 1 which is correct. 

For p >= 1 we look for terms where 2p-21-l = 1 or l=p-l, resulting in: 



C&O - (6-5) 
in case one prefers to express it the other way around. 
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It is easily seen from (|6.4f) that for large p: 



or: 
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lim=^ = l 



Renaat Van Malderen 

Address: Maxlaan 21, 2640 Mortsel, Belgium 

Author can be reached by email: hans.vanmalderen@bertholdtech. 



7 



